Antimatter production in supernova remnants 
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We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a super- 
nova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the 
trajectories of charged particles as a random walk with a prescribed diffusion coefficient, accelerating 
the particles at each shock crossing. Secondary production by CRs colliding with gas is included as 
a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The 
positron/electron ratio J-^+/J-^+^^- and the antiproton/proton ratio TpjTp^-p are a few percent 
and few x 10~^, respectively. Moreover, the obtained positron/electron ratio decreases with energy, 
while the antiproton/proton ratio rises at most by a factor of two above lOGeV. 

PACS numbers: 98.70.Sa, 95.30.Cq 
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I. INTRODUCTION 

Measurements of the antimatter fraction of cosmic rays 
(CR) provide not only insight into CR physics itself [l|, as 
e.g. their propagation in the galaxy, but are also valuable 
probes for cosmology and particle physics. In particular, 
the annihilation of dark matter (DM) leads to an equal 
injection rate of matter and antimatter particles into the 
Galaxy, while the CR flux from astrophysical sources is 
matter-dominated. A possible way to detect DM is there- 
fore to estimate carefully the expected antimatter fluxes 
from astrophysical sources and to search then for any 
excess 

The PAMELA collaboration presented recently results 
of their measurement of the positron fraction in CRs, 
which is rising rapidly from 10 to 100 GeV Q. At the 
same time, the antiproton ratio measured by PAMELA 
declines above 10 GeV Q, consistent with expectations. 
The conventional estimate for antimatter fluxes from as- 
trophysical sources uses as only production mechanism 
of antimatter the scattering of CRs on interstellar gas 
[l|. As discussed e.g. in the energy dependence of the 
Galactic diffusion coefficient, D (x with 6 = 0.5 — 0.6, 
is inconsistent with an increase of the antimatter fraction 
with energy. By contrast, the spectral shape of fragmen- 
tation functions leads quite naturally to such a rise in the 
case of DM annihilations or decays. 

The DM interpretation of the PAMELA excess faces 
however several difficulties 0: First, the required rate of 
positron production is larger than expected for a stable 
thermal relic. As a consequence, either the annihilation 
rate has to be enhanced by the dumpiness of DM or by 
non-perturbative effects operating at small velocities, or 
the DM particle should be unstable with the appropri- 
ate life-time. Second, in gauge boson or quark fragmen- 
tation, positron, antiproton and photon production are 
tied together and thus one has to postulate a DM parti- 
cle annihilating only into electrons and muons. More im- 
portantly, assuming antimatter production by diffusing 
CRs as the only astrophysical source for antimatter falls 



short: Since electrons lose fast energy, the high-energy 
part of the e~ + spectrum should be dominated by 
local sources as nearby pulsars, as pointed out already 
20 years ago Q. Moreover, electromagnetic pair cascades 
in pulsars result naturally in a large positron fraction to- 
gether with a "standard" antiproton flux. 

More recently, supernova remnants (SNR) were put 
forward as an alternative astrophysical explanation for a 
rising positron fraction Positrons created as secon- 
daries^ of hadronic interactions in the shock vicinity par- 
ticipate in the acceleration process and, according to Q, 
should thus have a flatter energy spectrum than primary 
electrons. It was estimated that the resulting positron 
fraction can explain the PAMELA excess and rise up to 
50% at higher energies 0, while subsequently a similar 
mechanism for antiprotons was suggested in Since 
shock acceleration in SNR is expected to be the main 
source for Galactic CRs Q, such a scenario has also im- 
portant consequences for the interpretations of CR data 
as, e.g., the boron-to-carbon ratio [l^- Additionally, sev- 
eral alternative explanations for the positron excess have 
been suggested: An inhomogeneous distribution of CR 
sources in the solar neighbourhood was put forward in 
Ref. 
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while the authors of Ref. [12| proposed an en- 
hanced secondary production in a cocoon-like region sur- 
rounding CR sources. 

The present work examines the production of sec- 
ondaryp and in SNRs, improving on previous stud- 
ies 0,11, [111 i'^ two respects: First, we use a Monte Carlo 
(MC) approach calculating the trajectory of each parti- 
cle individually in a random walk picture. This makes 
it easy to include interactions and the production of sec- 
ondaries. Second, our approach allows us to include the 
time (and spatial) dependence of relevant parameters de- 



^ Note that in the literature one often uses the notion of primary 
CRs in a more general sense for all particles produced by CR 
sources. Here we refer to the products of proton interactions in 
the source as secondaries. 
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scribing the evolution of a SNR as, e.g., the shock radius 
and its velocity, the magnetic field or the CR injection 
rate and to test their influence on the CR spectra. We 
should also stress what are not the aims of the present 
work: We do neither address the problem of acceleration 
from a microscopic point of view nor consider any feed- 
back of CRs on the shock or the magnetic field. Although 
the latter processes are important to obtain accurate CR 
escape fluxes, we shall show that our simplified treat- 
ment leads to an upper limit on the secondary fiuxes. 
Moreover we restrict ourselves to the production of an- 
tiprotons and positrons in the source, being the key issue 
in the proposal of Ref. 0, while we do not discuss sec- 
ondary production during CR propagation. 



II. SIMULATION PROCEDURE 

Shocks around SNRs are supposed to be collision- 
less, with charged particles scattering mainly on inho- 
mogeneities of the turbulent magnetic field. We model 
such trajectories by a random walk in three dimensions 
with step size lo{E) determined by an energy-dependent 
diffusion coefficient D. Diffusion close to the shock is 
usually assumed to proceed in the Bohm regime with the 
mean free path Iq proportional to the Larmor radius Rl ■ 
Thus D{E)(xE and 



D = 



= .fE 



3eB ' 



(1) 



where /b denotes the ratio of the energy density in the 
turbulent and in the total magnetic field. We neglect the 
coupling between CRs and the turbulent magnetic field, 
assuming that a layer with Bohm diffusion extends far 
enough into the up-stream region. For a constant mag- 
netic field, most low-energy CRs do not escape but are 
confined in the SNR, corresponding to an "age-limited" 
scenario for the CR flux from SNRs. 

We describe the evolution of the shock in the rest frame 
of the SNR. Then the (yet unshocked) up-stream region 
is at rest, ui = 0, and has the density of the surround- 
ing interstellar medium (ISM), pi = pisM- Assuming a 
strong shock with Mach number Ai ^ 1, the shocked 
down-stream region flows with the velocity V2 = 3vsh/R 
and has the density p2 = Rpi- Here, R denotes the com- 
pression ratio R = (7 -1-1)/ (7 — 1) = 4 for a mono-atomic 
gas with 7 = 5/3. We account for the flow, adding in 
the down-stream region on top of the random walk an 
ordered movement of the particle with velocity V2 that 
is directed radially outwards. Thus a particle trajectory 
evolves during the time step At ~ Iq/c as 



x(i + At) = x(t) + V2 At -dirsh - r) + Iq 



(2) 



where Ig denotes a random step, Vgh the radius of the 
spherical shock front at time t, and i?(a;) the step func- 
tion. 

Crossing the shock, particles are accelerated. We ne- 
glect that the relative energy gain ^ = {Ek+i — Ek)/Ek 



per cycle k depends on the angle of the trajectory to the 
shock front, and use for simplicity that on average for a 
non-relativistic shock ^ = ^(u2 — wi) = I'sh/c- For the 
position Tsh and the velocity Vs\x of the SNR shock we use 
the 7T, = case of the analytical solutions derived by 14 1. 
These solutions connect smoothly the ejecta-dominated 
phase with free expansion rsh oc t and the Sedov- Taylor 
stage oc t^/^. The acceleration of CRs is assumed to 
cease after the transition to the radiative phase at the 

time tmax- 

As the injected particles diffuse, electrons lose energy 
via synchrotron radiation, while protons can scatter on 
gas of the ISM producing secondaries that include an- 
tiprotons and positrons. Cross sections and the final 
state of pp-interactions are simulated using QGSJET-II 
p^ . while we use SIBYLL 2.1 [l^ for decays of unstable 
particles. 

The last ingredient for our simulation procedure is an 
injection model. To ease the comparison with the results 
of 01, we fix the electron/proton ratio K^p at injection to 
Kf.p = 7x 10"'^. As injection energy we use Eq = 10 GeV. 
In the first model used, the injection rate 



N ^rl^v'^^5{E~Eo)5{r~ny,) 



(3) 



is proportional to the volume swept out per tim e by the 
shock, i.e. a = 1 (thermal leakage model of Ref. [1^]). In 
the second model, the injection rate N is proportional to 
the CR pressure [201 ^'^'^ a = 3. In the case of model 2, 
the fraction of particles injected very early is significantly 
larger than in model 1. 

We use the following parameters to describe the SNR: 
We choose the injected mass as Mej = 4M0, the me- 
chanical explosion energy as Esnr = 5 x 10^^ erg, and the 
density of the ISM as tiism = 2cm^'^. The end of the 
Sedov- Taylor phase follows then as tmax = 13.000 yr [l3 |. 
For the magnetic field we use B = IpG and /s = 1, if 
not otherwise stated. 



III. NUMERICAL RESULTS 

In Fig. [1] we show the energy dependence of the pro- 
ton spectra in model 1 and 2. Additionally to the total 
spectra, the contribution of protons staying at tmax in 
the up-stream region is shown in red, while the spectra 
of protons advected to the down-stream are shown in 
blue. We do not show the spectra of electrons, since they 
have the same shape as the proton spectra apart from a 
somewhat lower cutoff energy due to synchrotron losses. 
While the total energy spectra at low energies agree well 
with a l/E"^ power-law, changing the injection model 
leads to large differences at high energies. The strong 
dependence of i?max on the injection model is expected, 
since the maximal energy is sensitive to how many par- 
ticles are injected early, when shock acceleration is most 
effective. However, one may wonder why a bump close to 
^-max exists in model 2, while model 1 seems to be well 
described by an exponential cutoff. 
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model 1 




E(eV) 

FIG. 1: Proton spectra (black) as function of energy for two 
different injection models and fsB — 1 /iG. Additionally the 
contribution of protons staying until tmax in the up-stream 
(red) and of protons in down-stream region (blue) are shown. 
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FIG. 2: Spectra of cosmic rays at tmax as function of energy 
E for different injection times t = 1, 300 and 3000 yr. 



Apart from energy losses and interactions that do not 
influenee the proton spectrum for the parameters chosen, 
deviations from a power-law can be introduced by 

an additional factor: The energy spectra of particles left 
behind the shock may be at any time t a scale-invariant 
power-law up to i? < £'max(i), while the total en- 
ergy of particles still participating in acceleration is con- 
centrated around £'max(i) as shown by the up-stream 
component in Fig. 1. If the latter carry a significant 
fraction of energy, as it happens in model 2, a bump in 
the total spectrum E'^dN/dE is visible. 

To demonstrate this effect, we show in Fig. [2] the spec- 
tra of CRs as function of energy E, obtained by injecting 
the same number of protons at different injection times 
i = 1, 300, and 3000 yr for model 2. Clearly, the peak 
becomes more pronounced with increasing i?max- The 
overdensity 5N of particles contained in the peak is small, 
but becomes pronounced in the E'^dN/dE plot since the 
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FIG. 3: The positron ratio J-^-ir / J-^+ (blue) and antipro- 
ton ratio iFp/Tp+p (red) in model 1 (dotted) and 2 (solid). 



associated energy 5E = ESN is significant. If this pic- 
ture is correct, then the CR spectrum of model 1 should 
also develop a peak close to -Emax at sufficiently late time 
t 3> tmax- We confirmed numerically that this is indeed 
the case. 

We switch now to the discussion of the secondaries 
produced. In Fig.[3]we show their energy spectra for the 
two injection models considered. In addition, in Fig. 2] we 
show for injection model 2 and fsB = the antipro- 
ton flux split into a part produced in the acceleration zone 

(A) and a part produced in the inner part of the SNR 

(B) . More exactly, we define the contribution A as all the 
secondaries that crossed at least once the shock. This 
contribution increases fast, since the time tacc primary 
protons stay in the acceleration zone and can interact in- 
creases as tacc oc D{E) oc E. Moreover, the inelasticity, 
i.e. the energy fraction transferred to all antiprotons is 
practically constant, (zp) w 0.02, in the relevant energy 
range, Eo/{zp) > 10^^ eV. Therefore it does not influence 
the shape of the antiproton flux. At i?b ~ 2 x 10^^ eV, 
the increase of contribution A stops, the total flux retains 
its approximate E~^ slope and stays small in contrast to 
the result of 01 . 

How can we understand this behavior and the maximal 
value of J-'p/J-"p_|_p? We may assume in a gedankenexperi- 
ment that in each pp ^ p + X interaction the most ener- 
getic antiproton carries away all the energy, Ep w {z)Ep 
with (z) « 1. Then interactions just convert part of the 
p into a p flux. But since p and p diffuse and are accel- 
erated in the same way, the total p flux is not affected if 
the p or the parent proton is accelerated. Hence the to- 
tal flux of antiprotons produced in the acceleration zone 
and inside the SNR, i.e. the sum of A and B, should be 
simply the proton flux scaled down by a constant fac- 
tor. In particular, the secondary flux of the species i 
is bounded by the proton interaction depth r and the 
(spectrally averaged) energy fraction (zj) transferred to 
i. The maximal conversion rate during the life-time of a 
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FIG. 4: The total flux of antiprotons together with the con- 
tribution A and B in model 2 as function of energy. 
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at 100 GeV given by r = cimaxi?niSM cr^p''' « 3 x IQ-^. 
The mean energy fraction of antiprotons (plus antineu- 
trons) is {zp) « 0.02, so we may expect a maximal ratio 
of J^p/Tp+p ^ (zp) T ^ 6 X 10~^. The obtained TpjTpj^p 
ratio shown in Fig. 0] is indeed close to this estimate. 

The same discussion applies to the case of positrons, 
with the sole exception that the primary electron flux is 
scaled down by the factor K^^ and that the energy frac- 
tion transferred to positrons is k, 0.05. The results 
of our simulation are shown in Fig. [3l confirming with 
the maximal value of T^+jT^+j^^^- < few % this sim- 
ple picture. Note that while above ~ 100 GeV the ratio 
J^e+ 1 ^e++e- from SNR starts to be larger than the con- 
ventional prediction using only secondary production on 
the ISM, it cannot explain the rise to J"e+ / J"e++e- « 10% 
at 10 GeV in the PAMELA data H. 



IV. DISCUSSION 

Both the absolute normalization of secondary positron 
and antiproton spectra and, more importantly, their 
spectral shapes obtained disagree with the results of 
0, d, [l3|. In the following, we discuss the reasons for 
these discrepancies. 



A. Secondary production 

One of the main advantages of the Monte Carlo ap- 
proach is the possibility to include interactions and the 
production of a multi-particle final state in an easy and 



self-consistent way. In contrast, the analytical calcula- 
tions of 0, S [31 involved certain approximations for 
the production of secondaries. In particular, a con- 
stant average energy fraction per single antiproton (or 



positron) 
and : 



6: = 
= 0.05. 



i)/{ni) was assumed, with ^p = 0.17 
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FIG. 5: The inelasticity (zp) and the multiplicity np of an- 
tiproton production as function of the primary proton energy 
E. 



In reality, it is the inelasticity, i.e. the energy fraction 



pp 



dE' 



„ E' d(jP^P{E,E') 



E 



dE' 



0.02 (4) 



transferred to all antiprotons, which is practically con- 
stant in the relevant energy range, Eo/{zp) > 10^^ eV. 
This behavior is illustrated in Fig. [5] where (zp) calcu- 
lated with our simulation using QGSJET-II is shown. 
(A similar result is obtained with SIBYLL.) The same 
holds true for the inelasticity of positron production, with 
(ze+) w 0.05. However, the multiplicity of secondaries 
rises relatively fast at those energies, as can be seen in 
Fig. El As a consequence, the average energy fraction per 
single antiproton (or positron) decreases strongly with 
energy. Therefore, assuming a constant average energy 
fraction per sin gle secondary, £^p = 0.17 and = 0.05, 
as done in 0, Isl. Il3l| , leads to an overestimate of by an 
order of magnitude for Ep ~ TeV. Moreover, the incorrect 
use of constant results in a wrong energy dependence 
of the secondary spectra. 



J 
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B. Comparison with the stationary picture 



Let us next discuss the differences between our results and the ones obtained in Refs. 0, H, [ll] that were based 
on a stationary picture and a simplified kinetic equation. In a stationary approach, it is more convenient to use the 
shock frame as reference frame, with the shock position at a; = and Ui as the flow speed upstream. The phase space 
density f{x,p) of CR protons is given by a power-law spectrum with exponent /3 — 3R/{R — 1), 



■P) 



pxm/D{p)^ a; < 
1, a; > 



Here, we have allowed for a spectral cutoff in the primary proton spectrum at p — Pn 
CRs is NcRiE,x) = 4Trp'^fcB.ix,p). 

Wc solve the transport equation for secondary cosmic ray phase space density 



(5) 



The energy distribution of 



dfsix,p) 
dx 

qs{x,p) 



dfs{x,p) 
dx 



Idu dfs{x,p) 



dEo NcrXEq, X) — nisM c , 

aE 



(6) 
(7) 



following the same steps as [T^ but avoiding any additional approximations for the source term qsix,p) for secondary 
particles. 

The ratio of the secondary CR flux Js{E) produced in the downstream region to the primary flux JcR(i?) follows 
then using a similar notation as Q as 



JcniE) 

where A{E) and B{E) arc however given by 



nisMc[A{E) + B{E)]^ 



B{E) 
A{E) 



ulE 



dz 



zf^'^daPP^''{Eo,z) 



E/E„ 
E 



dz 







dEoa';^;\E„)I{E„,l) 



Eo=E/z 



dE„a';f{Eo)I[Eo,^ 



E 
Eo 



(8) 

(9) 
(10) 



with 



dz . . , 




1 z^-^daPP^^'iEo, 



dz 



(11) 



In the left panel of Fig. |6l we show the spectral ra- 
tio Js{E) / JcB.{E) for the case of secondary antipro- 
tons as well as the partial contributions cnisMA{E) and 
cn'isMB{E). Wc employed in these calculations the 
same hadronic interaction model and the same source 
parameters {R = 4, nigM = 2cm~'^, ^max = IS.OOOyr, 
/bB^I /iG) as in our MC simulations, used a constant 
speed ui = 10* cm/s and varied the high energy cutoff 
^'max of the primary CR spectrum. The spectral index of 
the phase space density / of primary CRs has been flxed 
to (3 = A. 

First, we observe a much stronger dependence of the 
component A on the choice of E^^x as |8|. Inspection of 
Eqs. ([TU]) and pT|) shows that the dominant contribution 
to A{E) for E <^ i?max comes from the second integral on 
the r.h.s. of Eq. (ITUl) : for small e I{Eo,e) ~ e, and the 



main contribution comes from the l/z term inside the 
square brackets of Eq. pT|) . As a result, assuming the 
scaling picture^, A{E) diverges logarithmically, A{E) ^ 
D{E) \niE^,^/E), for E^,^ ^ oo. 

The authors of Refs. 0, replaced the factor l/z inside 
the integral of Eq. (|lip by its "characteristic value" 1 in 
front of the integral. Such a procedure would be permit- 
ted if l/z was replaced by its correct (energy-dependent) 
average with respect to the integrand of Eq. (|lip . The 
same procedure of "hiding" this logarithmic divergence 



^ In reality, Fcynman scaling is broken by the logarithmic rise of 
<j'j^p^{Eo) and by a small power- law like rise with energy of the 
central rapidity density of secondary particles. 
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FIG. 6: The antiproton ratio Tp/Tp+p (solid lines) for different i?niax (as indicated in the plot) as calculated for a fixed 
ui = fO* cm/s (left), and for ui derived via the relation D{Eraax) /u\ = tniax/20 (right) for the corresponding -Emax- Partial 
contributions to the ratio of the components A and B are shown as dashed and dotted lines, respectively. 



was copied by [13j. The incorrect averaging of the fac- 
tor 1/ z together with the misconception that the energy 
transferred to a single antiproton or positron per inter- 
action is energy independent are two major flaws in the 
analysis of 0,lll. 

Using the correct expressions, Eqs. ()10mip . the spec- 
tral shape of secondary cosmic rays can be almost arbi- 
trarily modified in the stationary approach, if E'max is 
treated as a free parameter. In contrast, using the usual 
relation for iSmax [HI 

i^(£^max)/w? -<max/20, 

the relative normalization of the A and B components is 
fixed. In such the calculated antiproton to proton 

spectral ratio, as plotted in the right panel of Fig. [51 
agrees qualitatively with our MC results and does not 
have the steep energy rise predicted in Refs. [tI [sI. [13|. 



C. Parameters determining the acceleration 
process 

In Sect, mil we have discussed only numerical results 
for constant JbB = 1/xG and one may wonder if a "bet- 
ter" choice of parameters can increase the antimatter 
fluxes. In particular, the analytical formula in the sta- 
tionary approach of 0, H, seem to imply that the 
contribution A of antiparticles produced inside the ac- 
celeration zone increases for weaker diffusion, i.e. larger 
D. However, the term D/vf regulating the importance of 
A limits also via tacc (x D/vf < tmax the maximal proton 
energy. Using a constant value Jb = 1/20 as in 0,H, 
thus reduces i?max by the same factor. 

The relative size of the partial contributions A and B in 
Fig-Hlcan be understood considering the relation between 
the time t^cc spent by protons in A, their final energy E oc 
face and thus the interaction depth ta in A as function 



of energy, ta o: face oc E. In particular, it takes all the 
life-time tmax of the SNR to accelerate protons to the 
highest energies, cf. the up-stream component in Fig. 1. 
For the component B, the optical depth tb of the parent 
proton is tb oc (tmax — ^acc), which explains why the two 
components A and B sum up to an approximately flat 
spectrum. On the other hand, the relative normalization 
of the components A and B in the stationary approach 
of 0, H, [13 been imposed by hand: treating i?max 
as an external parameter and increasing it relative to 
its natural value given by iacc = imax one enlarges the 
relative contribution of A, as illustrated in Fig. [51 Indeed, 
for i?max = 100 TeV and the parameters used by 0, ll], 
positron and antiproton energies E ~ -Emax imply as size 
of the diffusion zone Ldiff{E) ^ D{E)/v2 — 2kpc and as 
acceleration time t^cdE) ^ D{E)/v2 ~ 10^ yr. 

It is noteworthy that in the treatment of [l3[ the limi- 
tation due to the finite size of a SNR, Ldm{E) < L™^'^ = 
V2Tuia.x7 has been imposed, which significantly reduced 
both the relative normalization and the steepness of the 
energy-rise of the component A compared to the origi- 
nal treatment of 0]. Namely, L^iff{E) — L'^{^ was used 
for E > inbreak, with inbreak defined by the condition 
-^diff (-Bbrcak) = -^diff^' Howcver, a more severe constraint 
on the contribution of the acceleration zone to secondary 
antiparticlc spectra comes from the finite life-time of a 
SNR, which was not accounted for fully in the analytic 
treatments of 0,11, [11. 

One may try to justify the approach of 0, as a 

method to account in an effective way for an amplification 
and damping of the magnetic field. Namely, one can as- 
sume that the spectra of primary protons arc pre-formed 
in the early phase, when magnetic fields are strongly am- 
plified by non- linear effects [13, [13] , while the production 
(and re-acceleration) of secondary antiparticles is dom- 
inated by the contribution of the Sedov- Taylor phase, 
when magnetic fields are damped. 

We test this suggestion by considering a simple toy 
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coefficient drops above a certain energy and/or outside a 
sufficiently small radius Vsh+Sr to a value close to the one 
typical for the Galaxy, then CRs can escape up-stream 
instead of being confined down-stream. Clearly, this ef- 
fect reduces the contribution B. On the other hand, the 
bounds Te+/J-e++e- ^ few % and Tp/Tp+p < 6 x 10~'^ 
will become stronger, since also the time for interactions 
in the acceleration zone will be shortened. Since our max- 
imal values of J-p/J-p^p and J-^+ /J-'e++e- depend only on 
^max, which is lower in escape-limited models than in age- 
limited ones, we conclude that the contribution of SNR 
to the observed antimatter in CRs does not lead to pro- 
nounced rise of antimatter fractions and is smaller than 
estimated in earlier works. 



FIG. 7: Spectra of cosmic ray protons and positrons (scaled 
up by a factor 1000) together with the partial contributions 
A and B in model 2 for a time-dependent diffusion coefficient. 



model for a time-dependent magnetic field: assuming 
it to be strongly amplified in the early phase, with 
fsB = 100 before the transition to the Sedov- Taylor 
phase at = 240 yr, and using fsB ^ 1/20 /iC at i > t*. 
In Fig. [3 we show for this case the proton and positron 
spectra using the injection model 2. Protons that were in- 
jected early are accelerated up to few x 10^^ eV, while the 
bulk of CRs injected when the turbulent magnetic field is 
damped has a cutoff around 10^^ eV. The contribution A 
to the positron flux saturates at £^ ^ few x 10^^ eV, i.e. 
at the energy expected for /bB = 1/20 /iC In contrast 
to Fig. 131 the component B dominates now the high- 
energy end of the positron flux. It is easy to see the 
physical picture behind these results: the damped mag- 
netic field is unable to retain pre-accelerated protons of 
energies E > 10^^ eV in the vicinity of the shock front, 
which thus escape from the acceleration zone during the 
transition JbB/iiG = 100 — > 1/20 and escape far up- 
stream or are advected downstream. As a consequence, 
the re-acceleration of secondary antiparticles of energies 
E > 10^^ eV is only possible in the very beginning of 
the second phase, while at later times the antiparticles 
produced by CRs escaping up-stream and downstream 
are not longer accelerated. Although different and more 
realistic scenario for the dumping of the magnetic turbu- 
lence can be considered, it is obvious that the qualitative 
behavior of the spectra of secondary cosmic rays will not 
be significantly modified: A slower damping would pro- 
long the rise of the component A to higher energies while 
reducing its overall normalization (oc D/vf). It is worth 
stressing that the corresponding picture is essentially a 
non-stationary one and therefore the complicated inter- 
play of particle escape and re-acceleration can not be 
described properly in a stationary approach. 

Finally, we stress that the splitting between contribu- 
tion A and B is artificial and depends as well as the total 
flux on the definition of the escape fiux: If the diffusion 



V. CONCLUSIONS 

We calculated the energy spectra of CRs and their sec- 
ondaries produced in a supernova remnant using a simple 
random walk picture. In contrast to a previous predic- 
tion that the positron fraction J-^+ /J^e+_|_e- can rise up 
to 40%-50% for Kep = 7x IQ-^, we found that the ratio 
levels off at a few percent. This value corresponds to the 
expectation combining the interaction depth t « 3 x lO^'^ 
of a proton during the life-time of a SNR with the en- 
ergy fraction z ^ 0.05 transferred to positrons. Simi- 
larly, the antiproton ratio J-p /J'p+p does not rise beyond 
fcwxlO"^. Our results suggest that antimatter produc- 
tion in SNRs cannot explain the rise of the positron frac- 
tion observed by PAMELA. Since a rising antiproton 
fraction is neither expected from CR interactions with 
the ISM nor from pulsars, measuring J^p/Tp+p ^ 10"''' 
would be therefore a reliable signature for dark matter. 

The reason for the discrepancy with earlier works can 
be summarized as follows: Our MC results for the case 
of a constant magnetic field agree qualitatively with the 
stationary treatment, as one can see in the right panel 
of Fig. 6. In both treatments, one does not observe any 
pronounced energy rise of the fraction of secondary cos- 
mic rays produced by the sources. On the other hand, 
the left panel of Fig. 6 illustrates the main reason for our 
differences with Ref. Q: Treating the maximal energy 
^'max as a free parameter, one is able to modify arbitrar- 
ily the high energy behavior of the secondary spectra. 
Additionally, the approximation in the solution of the dif- 
fusion equation in Ref. Q infiuences the behavior of the 
spectra at lower energies, making them approximately 
independent of i?max ~ in contrast to the exact solution. 
Other differences between our MC treatment and the an- 
alytic one in Refs. 0, H, E^l have a much weaker impact 
on the obtained results. 

While nuclear fragmentation can be also treated in 
QGSJET, we have not included yet the nuclear decay 
chains required to predict, e.g., the boron-carbon ratio. 
Such ratios will provide after their measurement by the 
AMS collaboration a tool to clarify if the scenario pro- 
posed in Ref. (3] or the one discussed here is in better 
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